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We study one-dimensional quantum walks in a homogenous electric field. The field is given by a 
phase which depends linearly on position and is applied after each step. The long time propagation 
properties of this system, such as revivals, ballistic expansion and Anderson localization, depend 
very sensitively on the value of the electric field $, e.g., on whether $/(27r) is rational or irrational. 
We relate these properties to the continued fraction expansion of the field. When the field is given 
only with finite accuracy, the beginning of the expansion allows analogous conclusions about the 
behavior on finite time scales. 
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INTRODUCTION 

In this paper we analyze the surprisingly rich long time 
behavior of a simple quantum lattice system. Single neu- 
tral atoms in an optical lattice undergo a discrete time 
quantum walk by repeating the steps of first a unitary 
"coin" operation on the internal spin state, then a coher- 
ent spin-dependent transport by one lattice site and, fi- 
nally, the application of a phase depending linearly on the 
lattice site. The last step, which corresponds to a homo- 
geneous electric field, has recently been implemented by 
an acceleration of the lattice [1] . The crucial parameter 
is the phase difference per lattice site $. Depending on 

quite different long time behaviors have been observed 
[T], namely ballistic expansion, revivals, and Anderson lo- 
calization. The first challenge is hence to identify those 
$ for which different infinite time behaviors and related 
spectral properties occur. As we will see, this behavior is 
extremely discontinuous. It depends on whether <i>/(27r) 
is rational or not, which can never be controlled or ver- 
ified under finite error bars, and can therefore never be 
used to explain the differences seen in an experiment. 
However, if one does the experiment, and compares the 
behavior for <I'/(27r) a simple fraction and, say, the (ir- 
rational) golden ratio one sees quite clearly already on 
small time scales the qualitative differences predicted for 
the infinite time behavior (cp. Fig. [l] and Fig. |2] below). 
In this sense the rational vs. irrational distinction seems 
practically relevant after all. The second challenge ad- 
dressed in this paper is to resolve this apparent paradox, 
and to show how to extract qualitative information for 
large but not infinite times from a $ given with finite 
accuracy. 

The most famous case in which the dependence of a 
spectrum on the rationality of a parameter was shown is 
the Hofstadter butterfly [2^, which pertains to a two di- 
mensional lattice system in a magnetic field. The crucial 
property is the rationality of the magnetic flux through 
one unit cell. This system is analyzed in terms of the 
one dimensional Mathieu equation, whose variants have 
recently been analyzed with great success in terms of cer- 



tain dynamical systems [31. Interesting analogies to the 
present system are also provided by the quantum kicked 
rotator jj. Exactly the same system has been studied 
mostly in the rational case, namely numerically in [3 |5] 
and by exact diagonalization on a flnite ring [7j. In the ir- 
rational case Anderson localization [BHTU] is conjectured 
in [T. We agree that this is possibly the typical case as- 
suming $ to be a uniformly distributed random number, 
but we also exhibit irrational cases with a quite differ- 
ent behavior. Most extremely, we show the possibility of 
"hierarchical motion", which is characterized by an infl- 
nite hierarchy of time scales in which ever more perfect 
revivals alternate with ever larger ballistic excursions. 

An overview of our answers to the first challenge, the 
infinite time behavior, is given in Table |T] Some of the 
terminology will be explained in detail in the respective 
sections. The quasi-energy spectrum refers to the spec- 
tral properties of the unitary walk operator in the sense 
of, e.g., fTTJ Ch. VII. 2]. The properties in the table de- 
pend on the exact value of $ or, equivalently on its entire 
continued fraction representation. For an experimentally 
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TABLE I. Overview: The connection between the proper- 
ties of $, its continued fraction coefficients Ci, propagation, 
and the spectrum of the unitary evolution operator. Detailed 
explanations are found in the subsections for the respective 
columns. 



given field with finite accuracy, only an initial segment 
of this expansion can be reliably determined. Depend- 
ing on how one continues the expansion, any of the three 
classes is compatible with such data. However, the be- 
havior on a finite time scale can be reliable predicted 
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from the data. As in the infinite case the important in- 
put is the size of the continued fraction terms Ci, which, 
in particular, determine the quality of revivals. These 
relationships constitute the answer we can give to the 
second challenge. 



DEFINITION OF THE SYSTEM 

For definiteness and to fix notations let us explicitly 
define the class of systems under investigation. We con- 
sider particles on the ID lattice Z with a two dimensional 
space of internal degrees of freedom. Basis vectors are 
thus \x,a) with x € Z and a = ±1. The state depen- 
dent shift is defined by S\x, a) — \x + a,a). We consider, 
moreover, a coin operation C\x, a) — Cq^|x, /3) with 
a fixed matrix C € SU(2), i.e.. 



C 




(1) 



with |ap -f- l&p = 1. A standard case is the Hadamard 
walk with a = b = The third element will be a 

phase shift depending linearly on x, given by the opera- 
tor exp(ii$)|a;, a) = exjp(ix^)\x,a). The constant $ is 
referred to as the electric field. The overall time step is 
then given by 



cs . 



(2) 



The systematic way to generate such approximations 
is the so-called continued fraction expansion [T2|. For a 
number a; = > it is computed by taking the integer 
part cq — [xo\, setting xi = l/{x — Cq), and repeat- 
ing this process for xi to get ci, etc. The number x is 
then uniquely characterized by the sequence (co,ci, . . .). 
One gets an explicit sequence of approximating rationals 
ni/rrii, where the denominators TOj and numerators rii 
both satisfy the recursion — CiVi^i + ri_2 with ini- 
tial values hq = cq, n„i — 1, uiq = 1, m_i = 0. This 
is equivalent to an iteration in which at step i the de- 
nominator Ci is replaced by q + l/ci-)_i, leading to the 
typographical nightmare from which the continued frac- 
tions derive their name. The desired approximation of x 
of order 1 /m? is stated as 



mi 



< 



1 



Ci+im^ 



(4) 



RATIONAL APPROXIMATIONS 



Therefore, when $' ~ 2t: ni/rrii is a continued fraction 
approximation of $, the error in ([3| is of order (t/mi)^, 
so the i*^ approximation is valid roughly on the time scale 
of the denominator , and correspondingly longer if the 
next term q+i is large. 

Intuitively, equation ^ says that large terms Ci cor- 
respond to especially good approximations. For exam- 
ple, the continued fraction expansion of vr begins with 
(3,7,15,...) and the second approximant tt w 22/7 has 
an error < 1/(15 * 49) ~ 10""^. In this sense the "most 
irrational number" is the golden ratio ip = {1 + \/5)/2 
with expansion (1, 1,1,.. .). 



We will analyze the behavior of Wg, and its iterations by 
comparison with similar walks W^', where $' = 2Tm/m 
is a rational approximation to $. The basic tool is an 
estimate for the difference between two states evolved 
from the same initial state by walks with slightly dif- 
ferent fields $ and As initial state we choose a 
unit vector ijj which is non-zero only for lattice sites x 
with < L. Then for one step IjW^^V' — < 
max|a;|<L I exp(za;$)— exp < <I>'|. After i steps, 
adding the deviations, and observing that the localization 
region of -0 increases by 1 in every step we get: 

\\Wl^ - Wi,ij\\ <*-{t + 2L-l)\^- $'1 . (3) 

So, for example, if $' approximates <i> to within e = 10^^, 
and the particle starts at the origin {L — 1), we can 
use W^itp to predict the behavior of W^ip for up to 
£~i/2 « 100 steps. The square root makes decimal 
approximations (or approximations in any other digital 
base) quite useless. Instead, one can look for denomina- 
tors m such the approximation error is not of order m~^, 
as can always be trivially achieved, but of order m"^ or 
better. 



THE RATIONAL CASE: REVIVALS AND 
BALLISTIC EXPANSION 

Let us begin with the rational case, in which <1> = 
27rn/m, where it is understood that n,m do not have 
common factors. The phase factor 77 — exp(i$) is then a 
primitive m}^ root of unity, i.e., 77™ — 1, but this holds 
for no smaller power. There are two ways to reduce the 
analysis to the case of standard, translationally invariant 
walks |13ffT5] . One is to use the fact the field phase fac- 
tor is periodic with period m, and hence always group m 
sites together into a "supercell" of "coin dimension" 2™. 
This spatial regrouping is carried out in [T] , giving a band 
structure with a 2™-branched dispersion relation. Here 
we will use instead a temporal regrouping, i.e., consider 
the operator W^, which likewise commutes with trans- 
lations. This is a walk with two coin states, but allowing 
for steps over up to m sites. The key fact about this walk 
is that, even for moderately large m, it hardly moves, in 
the sense that the walk operator is close to the identity 
operator (up to a phase). Phrased as a statement about 
the ungrouped walk W^, we get a revival of the initial 
state after m steps. The precise statement is the follow- 
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ing revival theorem: 



Wl'^ + 1 1 



2\a\ 



odd 



op 



(5) 
(6) 



Hence for odd m any initial state is reproduced after 2m 
steps with accuracy (trace norm distance) 4101™, which is 
exponentially small in m, since \a\ < 1, except in trivial 
cases. For the Hadamard walk we get an error 2^'"/'^. 
The different behavior of the even and odd cases is to 
be expected for coined walks, since the amplitude for no 
jump or a jump by an even distance is zero. Hence the 
revival of a state localized at a point is only possible after 
an even number of steps. 

The proof of the revival theorem will be given in the 
appendix, and could also be based on the eigenvalues 
given, without proof, in |7]. It yields also the dispersion 
relations for concretely cosa;-|-(fc) = 



a™ cos{mk) 
-cos(mfc) + 



odd 



(7) 



The upshot is that the spectrum is absolutely continuous, 
and transport is ultimately ballistic. However, because 
of the revivals it may take a very long time for the bal- 
listic regime to be reached. This is shown in Fig. [T] for 
walks starting at the origin. In each case we show two 
aspects as a function of time, the root mean square devi- 
ation of position a{t) = (.t^)^/^ and the probability p{t) 
of return to the origin. Perfect revivals are identified by 
the conditions a{t) = 0, resp. p{t) = 1. The two fields 
chosen in the figure are $i = 27r/5; for which ^ pre- 
dicts a fairly weak revival at t = 10 (p(10) > .64), which 
rapidly goes over into ballistic expansion. The envelope 
ofp(t) is well approximated by a Bessel function The 
denominator <I>2 = 27r x 51/256 for the second field is 
much larger, so the revival predicted by ^ is exponen- 
tially sharper: p(256) > 1 — 10^^^. The full evolution up 
to t = 256 therefore repeats roughly 10^^ times, i.e., un- 
til the revival errors accumulate sufficiently, to make way 
for ballistic expansion. This would be true for all fields of 
the form $ = 27m/256, independently of the numerator 
n. However, the evolution up to t = 256 does depend on 
the numerator [5]. We can get some information about it 
from a "rational approximation of the rational" 51/256: 
Its continued fraction sequence is (co,ci,C2) = (0,5,51) 
and 1/5 — $i/(27r). Therefore, we know that the initial 
segments of these evolutions coincide. More precisely, 
from (§ we get \\W^^^P - W^^^I^W < t{t + l)/tl with 
to « 20. The comparison for the displayed quantities 
is even more favorable, as shown in the insert of Fig. [l] 



THE ALMOST RATIONAL CASE 

The same idea, i.e., combining (|3|,(|4]), can be used with 
the revival theorem ([s]) resp. ([6| to get revival estimates 




FIG. 1. Revivals for two walks with rational fields. Top 4>i = 
27r/5 and bottom <1>2 = Stt x 51/256. In both cases the right 
panel shows the root mean square of the position, and the left 
panel the probability p{t) to be back at the origin at time step 
t. Since $1 — $2 ~ 5 * 10"^ the initial segments {t < 100) of 
the top and bottom graphs almost coincide (see black frames 
and insert). The revivals predicted by the revival theorem are 
at t = 10 and t = 256, respectively. 



for arbitrary irrational $: 



M^4"'V + (-l)"'/V 



< 



in 

Ci+l 
TT 



o 



o 



L 



L 



odd 



m, even. 



(8) 

That is, if Q^+i 00 along some subsequence ia we get 
an infinite sequence of sharper and sharper revivals at the 
corresponding times m^^ or 2mi^ . From each of these re- 
vivals we also get a more precise repetition of the entire 
history up to that point. This shows that the walk opera- 
tor now has no absolutely continuous spectrum, since for 
initial vectors in the absolutely continuous subspace the 
expectation for any finite dimensional projection goes to 
zero by the Riemann-Lebesgue lemma. 

When the Ci grows sufficiently fast, we get hierarchical 
motion on a hierarchy of time scales m^, each of which 
is associated with one approximation rii/mi. Sharp re- 
vivals, up to errors £i with e^— >0, alternate with large 
excursions, whose length is given by a sequence of in- 
creasing intervals 7^. The walk will have the property 
that for any state initially supported by li there is a 
time when the probability to be back in li is less than e^. 
The sequences Si and li can be prescribed arbitrarily. 

We now show how a walk with these properties can be 
constructed by judiciously fixing the sequence (cq, ci, . . .) 
step by step. Suppose the sequence is constructed up to 
Ci- Then a revival is predicted by (|8]) at or 2mj. 
Clearly, by choosing Ci^i sufficiently large we can make 
the error in ([sl less than e^. Now the rational walk with 
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$i = 27mi/mi will eventually turn ballistic. That is. 
there is an "excursion time" at which any initial state 
supported on li will be back in li with probability less 
than £i/2. Choosing c^+i even larger if necessary we 
can ensure by ([3| and ^ that the approximation will 
hold up to this excursion time with error less than ei/2. 
Hence the constructed walk will have the desired excur- 
sion property. 

We remark that in this construction the will have to 
grow very fast indeed, since after a sharp revival the bal- 
listic stage is only reached after many repetitions. The 
fields <I> arising from our construction are therefore not 
typical (the set has Lebesgue measure zero). Hierarchi- 
cal motion implies that the spectrum of W$ is purely 
singular continuous. In fact, we have already excluded 
absolutely continuous spectrum as inconsistent with in- 
finitely many revivals. On the other hand there can be 
no eigenvalue (point spectrum). Indeed an eigenvector 
would be essentially (up to some e) supported by some 
li- Since an eigenvector remains constant in time under 
the walk this contradicts the excursion property. 



THE VERY IRRATIONAL CASE 

Numbers with bounded continued fraction sequence Ci 
do not have exceptionally good rational approximations. 
For example, for the golden ratio ip with Ci — 1 for all i, 
the bound (|8| makes no revival predictions at all. Fig. [2] 
shows the behavior of root mean square position and re- 
turn probability for = 2Trip. It shows many peaks, not 
just at the continued fraction denominators, and no ten- 
dency for transport or decaying return probability. This 
is precisely the behavior one would expect from a quan- 
tum system with pure point spectrum. Indeed a numer- 
ical study readily confirms this guess. 




FIG. 2. Root mean square of position and return probability 
for the golden ratio field <1> = 7r(\/5 — 1). There are many 
returns, and the initial state has a substantial overlap with a 
bound state, so the return probability is bounded away from 
zero. Vertical lines: Times appearing in the revival estimate, 
even though the estimate is trivial. 

More precisely, we take the system with periodic 
boundary conditions, diagonalize it, and pick the eigen- 
function with the smallest (x^)- This eigenfunction de- 
cays exponentially (as long as one does not approach the 



boundary), and converges as the ring size is increased. 
The same function is obtained by iterating the trans- 
fer matrices [9] for the solution of the eigenvalue equa- 
tion W^ip = exp(iijj)ip. The quasi-energy uj is chosen as 
w = $/2 as this ensures the symmetry of the eigenfunc- 
tions around the origin. Therefore, if we start with an 
arbitrary vector at some very negative site x — —N , the 
iteration will naturally pick the expanding branch, but 
on passing x — Q will decrease again to reach exactly 
the starting value at a; = -l-A^. When N is large, this 
computation has to be carried out with high numerical 
precision (see Fig. [s]). Setting the function to be zero 
for |a;| > N yields an approximate eigenfunction, which 
solves the eigenvalue equation with a precision (e.g., 2400 
digits for the case shown in Fig[3|, which in ordinary nu- 
merical work would be considered extraordinary high. 
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FIG. 3. Position probability distribution for the eigenfunction 
of the walk with $ = 2Tiip. Left: Unsealed plot. Right: Loga- 
rithmic plot over the range [— A^, A'^] for which the procedure 
described in the text was applied. Floating point accuracy for 
this computation was set to 8000 decimal digits. The slope 
of the right graph is the inverse Anderson localization length 
^ .301. 



Knowing one exact eigenfunction is knowing a com- 
plete set, since we can generate the required number 
per site by shifts and staggered sign changes ij{x) i— >■ 
(— l)^-!/)(a;). We emphasize that in spite of the high pre- 
cision the numerical evaluation does not prove the eigen- 
value equation, since for larger N the behavior could be 
different, and indeed will be different for rationals close 
enough to <I>. Hence we have shown Anderson localization 
"only for all practical purposes", meaning for runs with 
less than ^ g^gpg^ which for nanosecond steps will 

be incomparably larger than the age of the universe. 

It is paradoxical that proving localization for a con- 
cretely specified irrational is very hard, whereas it is pos- 
sibly much easier to show that localization holds for all 
$ outside a set of measure zero. There are some indica- 
tions for the latter statement to be true. Indeed one can 
repeat the above eigenvector computation for random <I> 
and finds that the inverse localization lengths are always 
around .301 with a variance that decreases when N and 
the numerical accuracy are increased. In contrast, the 
small scale behavior (left panel in Fig. |3| of the eigen- 
function does depend visibly on $. We note that this 
conjecture does not contradict the markedly different be- 
havior for rational and almost rational $. since these sets 
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are of measure zero. 



SUMMARY AND OUTLOOK 

We have established three typical kinds of long time 
behavior of electric walks as a function of the field 
This dependence is extremely discontinuous in the sense 
that an arbitrarily small change in $ may change the 
type of propagation behavior in the infinite time limit. 
Nevertheless, on a given time scale one can make precise 
predictions in terms of an initial segment of the continued 
fraction approximation of $. Anderson localization holds 
for very irrational numbers, and we conjecture that it also 
holds for random numbers $ with probability one with a 
non-random localization length. 



APPENDIX 

The revival theorem is stated in equations ^ and ([6| 
in the body of the paper. We repeat it here for reference: 

Theorem 1. Let $ = 2tt— such that n and m are rel- 

m 

atively prime integers. Then the electric walk operator 
satisfies the revival relations 



\\wl^ 



I =2\a\ 

W op ' ' 



= 2|ar'/2 



m odd (9) 
m even . (10) 



In order to prove this theorem, we analyze the opera- 
tor W$ in momentum space. That is, we consider wave 
functions on the Brillouin zone with values in the coin 
space C^. Operators commuting with translations are 
then precisely those which act by multiplication with a 
matrix depending on the quasi-momentum k. In partic- 
ular, the coined walk in zero field acts by Wii{k)'ip{k), 



Wo{k) = S{k)C : 



(11) 



where S{k) = expika^. The operators e"* act in mo- 
mentum space as a shifts of the quasi momentum, i.e., 
(e'^*'0)(fc) — il){k + where the addition is understood 
mod27r, as usual for quasi momentum, so in particu- 
lar ((e"''')'"V')(^) = Applying W$ therefore gives 
{W^'tj^){k) = Wo{k + <^)^{k + $). With m applications 
the momentum is shifted by to$, which by assumption 
is a multiple of 27r, and can hence be left out. This 
precisely reflects the anticipated translation invariance 
of WJ^ . Hence, this operator is represented by the k- 
dependent 2 x 2-matrix 

W^'{k) = Wo{k + $) • Wo(fc + 2$) • • 
= S{<^)WQ{k)S{<^f ■ ■ ■ S{^Y 



■ W{k + m<I>) 
Wo(fc) . (12) 



By definition, the dispersion relations associated with 
such a walk are the functions u)± (k) such that the eigen- 
values of W^{k) are exp(ia;±(fc)). Since we are dealing 
with 2x2 matrices and det{W^ (k)) = 1 the dispersion 
relations depend only on the trace of W^{k), more pre- 
cisely, we have 



2cosw±(fc) = tr(iy^(A:)) . 



(13) 



Evaluating this trace could be expected to be a hard task 
involving a sum over 2™ oscillating terms. Surprisingly, 
however, there is a simple formula for this type of expres- 
sion. We state it slightly more generally than needed. 
The proof of the Theorem will be completed after the 
proof of this Lemma. 

Lemma 2. Let m G N, and rj he a primitive m*^ root of 
unity. Consider the matrices 



C ^ 




id R = 



(14) 
(15) 



V \ 

Of]-')' 

and set 

t™(C) = tr(Ci?"Ci?i • • • Ci?"-i) . 
Then for odd m: 

T™(C) = a™ + d™, 

and for even m: 

r„(C) = - (a™ + d™) - 2((-ad)"/2 _ det(C)"/2) 

Proof of the Lemma 



For any n (not necessarily n = m), consider the expres- 
sion T„(C) written out as the sum. 



I2l'i ■ ■ ■ ^Inll^iji 



(16) 



Then with each term we can associate a closed path ii 
^2 —>••.—). i,j of length n on the graph 



a 




This means that every summand contains the same num- 
ber of h and c factors, so t„(C) depends only on the 
product be. Hence we can write it as a polynomial in the 
variables a, d and 6c; 



Tn{C) = Y,^rsa'd' (be) 



n—r—s 



)/2 



(17) 
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with suitable coefficients 7^5, where the degree of (6c) 
follows from the homogeneity of degree n of Tn{C) in all 
four variables (a, b, c, d). 

Now assume that m is odd, set n = to in (17 1, and 
substitute C ^ RC . This does not change the trace 
because i?™ — 1 = and we can cyclically shift the 
factors CR^ under the trace. On the other hand, this 
substitution replaces a t-^ rja and d 7]^^d. Thus the 
r, s-term in (17 1 gets multiplied by 77'^''"''^. Since r,„(C) 



is unchanged, the uniqueness of the polynomial coeffi- 
cients tells us that the only contributing monomials can 
be those with r = s mod to. Since < r, s < to it follows 
that either r = m, s = 0, r = 0, s = m, or r — s. How- 
ever, the last possibility is ruled out, because the overall 
degree of the trace as a polynomial in (a, 6, c, d) is to, 
and hence odd. Therefore the only non-zero coefficients 
are 70m and 7om- Clearly, these coefficients are just the 
product of the respective powers of R, namely 77^ with 
N = jy^=o ^ ~ m{m — l)/2. Since (to — l)/2 is an inte- 
ger, we have 7]^ = 1, which entails the expression given 
in the Lemma. 

The argument in the even case is similar. Let us set 
TO = 2m'. Our first aim is to compute Tm'{C). We 
introduce the operators 



Xf — CR ■ ■ ■ CR 



'+m'-l 



(18) 



i.e. products involving a contiguous block of only half 
the terms in the trace. By definition. 



tr{X,) =T„ACR') 



(19) 



From 77™ — 1 = {rj™ )^ we get 77™ — ±1. But since 77 is a 
primitive to*'' root, 77™ = — 1, and R™ = — I. Inserting 
this relation into the trace of Xi we find the relation 

tr(Xi) = tr(Ci?Ci?2 • • • CR"''-^CR"'') 

= tr{CR"''CRCR^ ■ ■ ■ CR"''-^) 

= - tr(Xo) . (20) 



T„ACR) + T„ACR^)^0 



(21) 



With (17 1, we see that the polynomial on the left hand 
side can only vanish if, for all r, s. 



1=0. 



(22) 



Since —to' < r — s < to', and the only integer powers t e 
[—to', to'] with 77* = —1 are t = ±to', we conclude that, 
as in the odd case, 7^^ 7^ only for (r, s) = (m',0) or 
(r, s) = (0,to'). Also as in the odd case 7m'o = rj^ with 

N = Efclo^ ^ = "^'(^' - We know that 77"^' = -1 

so we can write 7„i'o = and 7oto' = i^jY" with 

j — ±i. Which sign is correct depends on the root 77: 
We can write the primitive root as 77 = exp(27rm/TO) = 



exp{iTTn/m') for some odd n coprime with to'. Then 
77^ = exp((i7r77/2)(TO' — 1)), i.e., j = i". To summarize: 



Tm'(C)=j a +[-.]) d 



(23) 



Later on we will need only the square of this trace, so the 
ambiguity about the sign of j drops out: 



Tra'iCf = (-l)™-i(a™-f d"^) -}-2a'" d 



(24) 



Now the full trace we want to evaluate is t„i(C) = 
tr(XoX„i') and since X^i = (—1)™ Xq this can be done 
by the identity tviXf - tr{X^) = 2dct(X) for 2x2- 
matrices: 

tr{XoXm') = i-ir' ty{X^) 

= (-l)"'((trXo)2-2dct(Xo)) 

= (-l)"'(r„,,(C)2-2det(C)™') 

= -(a™ + rf™) - 2((-arf)"' - det(C)'"') 
which finishes the proof of the Lemma. 

The dispersion relation 

Coming back to the walk in a rational electric field 



we can evaluate (13 1 with il2« using the Lemma to get 



cosaj±(fc) — ^Tm{Wz{k)). Thus in the expressions of the 
Lemma we have to make the replacements a 1— exp(7fc)a, 
d I— exp(— ifc)a and det C 1. Setting a = \a\ exTp{iko), 
we get 



cos uj±{k) 



(25) 



laj™ cos(m(/c + fco)) to odd 

-|a|"cos(TO(fc + /co)) + (-l)"/2+^(l- |a|™) to even 



Proof of the revival theorem 



Intuitively, we will use that, according to the dispersion 



relations (251, all w±(fc) must be close to ±i in the odd 
case and close to (— 1)^to/2) in the even case. To get 
the precise norm expressions, consider a normal opera- 
tor A (i.e., = 0), which commutes with transla- 
tions, so its action on momentum space wave functions 
is {A7P){k) = A{k)ip{k). Then P|| = max^ ||A(fc)|l = 
maxk^± |Q!±(fc)|, where a±{k) are the eigenvalues of A{k). 
We can apply this to the operators in the revival theorem. 
In the odd case, 



wi 



I — max 

Hop fe,± 



max 

fc,± 



= max 2 cos uj± (k) 

k,± 

= 2\ar 



(26) 



7 



where at the last Hne we used that the maximum over k 
is attained at fc = — fco- Similarly, in the even case we 
use that \\A\\^ = \\A^A\\ and get 



— max 

op k.± 



+ (-1) 



m/2 



max 
fc,± 



(2 + 2(-l)"/2cosw±(fc) 
= (2(l-(l-|ar)) + 
+2|a|™ max(-l)"/2 cos(TO(fc + ko)) 

k 

= A\ar (27) 
This concludes the proof of the revival theorem. 
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